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I/-) ABSTRACT 

(~) Wide field near-infrared observations and Spitzer Space Telescope IRAC observations 

of the DR21/W75 star formation regions are presented. The photometric data are 
used to analyse the extinction, stellar content and clustering in the entire region by 
using standard methods. A young stellar population is identified all over the observed 
field, which is found to be distributed in embedded clusters that are surrounded by a 
distributed halo population extending over a larger projected area. The Spitzer/IRAC 
data are used to compute a spectral index value, a, for each YSO in the field. We use 
O i 1 these data to separate pure photospheres from disk excess sources. We find a small 

fraction of sources with a in excess of 2 to 3 (plus a handful with a ~4) , which is much 
higher than the values found in the low mass star forming region IC348 (a ^ 2). The 
sources with high values of a spatially coincide with the densest regions of the filaments 
, and also with the sites of massive star formation. Star formation is found to be occuring 

in long filaments stretching to few parsecs that are fragmented over a scale of ^lpc. 
The spatial distribution of young stars are found to be correlated with the filamentary 
nebulae that are prominently revealed by 8/im and 850 /im observations. Five filaments 
r jJJ ■ are identified that appear to converge on a center that includes the DR21/DR21(OH) 

regions. The morphological pattern of filaments and cluster i ng co mpare well with 
numerical simulations of star cluster formation bv lBate et alJ l)2003|k 
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1 INTRODUCTION 

Wide-field near-infrared (NIR) and mid-infrared (MIR) sur- 
veys are potentially very useful for mapping the distribution 
of young stars in Giant Molecular Clouds (GMCs). In partic- 
ular, the launch of the Spitzer Space Telescope has brought 
about a revolution in such studies, and several nearby, low- 
mass star forming regions have been ana lysed recently (ex 
Ijorgensen et alJl2006l : lHarvev et alj|200(f) . With the advent 
of wide-field imagers on mid-sized, ground based telescopes 
like the U.K. Infrared Telescope (UKIRT), the deeper near 
infrared data can be effectively combined with the Spitzer 
Space telescope observations and the techniques employed 
on nearby low mass cores can then be applied to relatively 
distant (d~2-3 kpc), high- mass star forming regions. 

Two such regions are DR21 and W75N in the Cygnus X 
HII complex. Traditionally observed separately, DR21 and 
W75N are associated with massive, dense cores separated 



by about half a degree on the sky dWilson fc Mauersbergerl 
Il99d: IShirlev et all 12003ft . Both sources are associated 
with ultra-compact HI I regions (lHaschick et alJ Il98ll : 
ICveano wski et aT| |200 3ft. wa t er, OH and/or methanol 
masers jPlambeck fc Mentenl Il990t iHunter et alJ Il994l : 
iTorrelles et al.lll997t iKoean fc Slvshl Il998h. and extremely 
mass i ve bipolar outflows jMoore et al.lll99lt iGarden et alJ 
119911: iDavis fc SmithllT999 : IShepherd. Testi fc StarkH2003Tl . 
DR21(OH), a third site of massive star formation, 
lies ~3' north of DR21. Althou gh it too is associ- 
ated with intense maser activity jKogan fc Slvshl Il998l : 
iManeum. Wooten fc Mundvlll992li . unlike DR21, it is not 
detected at radio wavelengths and its bipolar outflow is 
not seen in H2 emission. All these regions are bright 
at Far-Infrared (FIR) wavelengths; DR21 is estimated to 
have an FIR lu minosity of 1.5xl0 5 Lp) and the W75 re - 
gion 5-8x10 4 Lq llHarvev et al 1ll977UCampbell et aljll982Tl . 
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DR21(0H) probably harbours multiple sources at a younger 
stage of evolution. DR21 and DR21(OH) coincide with just 
two peaks in a chain of molecular cores, traced in molecu- 
lar line and submillimeter co ntinuum emission, that extends 
north-south over at least 12' JChandler. Gear fc Chinll993t 
IVallee fc Fiegell200ft . lKetol fll990() found collapse signatures 
in two of the main protostellar condensations coinciding with 
the UCHII regions DR21 and DR21 D and estimated masses 
of 270Mq and 20Mq respectively. All of the above suggest 
that DR21/W75N is a region with intense massive star for- 
mation activity. 



The DR21/W75 region has recently been observed with 
Spitzer and the wide-field camera, WFCAM, at UKIRT. 
An overview of the Spitze r IRAC and MIPS data has 
been given bv | Marston et al] 1:2004); Smith ct al. (2005J) and 
IPersi. Tapia fc Smitbl l)200ff) examine the MIR photometry 
of the massive, embedded young stars associated with just 
the DR21 and W75N cores, respec tively. The WFCA M ob- 
servations have been presented bv lDavis et alJ fepQrj. here- 
after Paper I), where narrow-band images in H2 1-0S(1) 
emission have been used to search for outflows and jets 
across a 0.8°x0.8° field. In Paper I the Spitzer/IRAC pho- 
tometry is used to search for candidate outflow sources; an 
850 /jm mosaic is also presented, which traces the W75N 
cloud and the chain of very young (possibly pre-stellar) 
star-forming cores running through DR21 and DR21(OH). 
A colour image showing the field observed with Spitzer is 
shown here in Fig. 1, along with contours of the 850 /im dust 
continuum emission, for reference. The relationship between 
the many outflows, the embedded young stars, the UCHII 
regions, masers and molecular cores is discussed in some de- 
tail in Paper I. In this paper - Paper II - we step back and 
take a more global view of DR21/W75: the Spitzer/IRAC 
photometry and UKIRT/ WFCAM JHK data are used to 
characterize and map the distribution of young stars across 
the field. We derive extinction maps, surface-density maps, 
and plot colour-colour and colour-magnitude diagrams for 
some 43,000 WFCAM sources and 1580 Spitzer targets. We 
use these results to study the nature of the clustered and dis- 
tributed populations, and their relationship to dense cores, 
filaments, and the locations of massive stars. 



DR21 and W75N are associated with clouds with aver- 
age Local Standard of Rest (LSR) velocities of ~-3kms _1 
and ~+9kms~ 1 , respectively; a small amount of foreground 
gas at Vlsr ~ + 9kms~ may also lie in front of DR21. 
iDickel et all <1978h propose that the DR21 and W75N cores 
are interacting, based on their modest-resolution mapping 
and an observed velocity gradient across W75. However, in 
our higher-resolution 850 /im maps in Paper I we fail to de- 
tect dust continuum emission between W75N and the 12'- 
long DR21 filament. Moreover, a distance of ~ 2 kpc is usu- 
ally adopted for W75N, while 2-3 kpc is a s sumed for DR21 
and pR21(OH1 iCenzel fc Downesl Il977t ICampbell et all 
11981 iFischer et alJll985t lOdenwald fc Schwartzlll993h . The 
higher value gives a more realistic population of massive 
stars in our analysis below, so we adopt 3 kpc as a general 
distance to DR21, DR21(OH) and W75N regions in this pa- 
per. 



2 DATA SELECTION AND ANALYSIS 
2.1 JHK Data Selection 

WFCAM observations of a 0.8° x 0.8° field, known as a WF- 
CAM "tile", were obtained in Service Observing time at 
UKIRT in May and June of 2005. Data through Mauna Kea 
Consortium broad-band J, H and K filters, and through a 
narrow-band H2 1-0S(1) filter were secured; the H2 data are 
presented in Paper I, where the details of the WFCAM ob- 
servations and data reduction are also discussed. The WF- 
CAM field encompasses all of the region covered by the 
Spitzer data displayed in Fig. 1. 

Catalogues of WFCAM point source photometry in the 
J, H and K bands were provided by the WFCAM Sci- 
ence archive maintained by the Cambridge Astronomical 
Survey Unit (CASU). JHK data were retrieved from the 
merged JHK catalog. The following four constraints were 
used when querying the catalog: a) magnitude errors of less 
than 0.2 mag in each band, b) a merged class identified as a 
star in the catalog, c) a probability that the point is a star, 
given by the pStar parameter, of greater than 0.999, and d) 
coordinates matched to better than 0.4" between the JHK 
bands (x and r\ within ±0.2" in each band). Regions of star 
formation contain large numbers of embedded sources which 
may not be visible in the J band. Therefore, the unmerged 
catalogs were used to obtain data in the H and K bands sep- 
arately to find matching sources only in these two bands. 

The JHK merged catalog with the above constraints 
yielded 42792 point sources in the region; the HK merged 
catalog contained an extra 16883 resulting in a total of 59675 
sources detected in the HK bands alone. These two catalogs 
have been used for the photometric analysis presented in 
this paper. The WFCAM photometric system i s based on 
the s ystem of Mauna Kea Consortium Filters llDve et alJ 
2006). We have maintained the same system for plotting 
the data points, although we have transformed the data of 
the reference curves (which are in Bessel-Brett system) to 
that of the MKO system. 



2.2 Spitzer Photometry 

Archival Spitzer Space Telescope data were retrieved using 
the Leopard software for the analysis presented here. We 
obtained the Post-Basic Calibrated Data (PBCD) images 
from two programs, namely PID-623 (IRAC Campaign W 
(SV-3), PI: Giovanni Fazio) and PID-1021 (DR21 and its 
Molecular Outflow, PI: Anthony Marston). These programs 
together cover almost the entire WFCAM field of view pre- 
sented here and in Paper I, encompassing the DR21, W75N, 
the Diamond ring regions to the south west of DR21 (see 
Fig. 1), as well as new sites of star formation to the west of 
W75N, and to the west of DR21 (labeled L906E in Fig. 1 
and in Paper I). Data were obtained with both the MIPS 
and IRAC cameras, although here we discuss only data from 
the latter. IRAC images in four filters, known as band 1 
(band center wavelength ~3.6 /im), band 2 (4.5 /xm), band 3 
(5.8 /im) and band 4 (8.0 11m) were obtaine d. The IRAC in- 
strument is described bv lFazio et alJ i2004f) . Both programs 
used Full Array Readout with 2 sec exposure times. Mo- 
saics of the individual frames were made using STARLINK 
routines and aperture photometry was extracted using the 
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Figure 1. A colour image of the full DR21/W75 mosaic observed with the Spitzer Space Telescope, composed from images in 3.6 /xm(blue), 
4.5 /im(green) and 8.0^tm(red). Contours represent the 850 fim emission mapped using SCUBA on JCMT. The sketched lines are meant 
to indicate the large scale filamentary emission. 
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Figure 2. Near-infrared extinction map shown in grey scale. 
Overlaid contours represent nearest neighbour stellar density en- 
hancements as derived from the 4.5^tm IRAC image point source 
photometry data. The contour levels show star densities of ~ 30K, 
40K, 50K, 70K, 100K, 150K and 200K stars degree" 2 

APPHOT tasks in IRAF. Images and contour plots in three 
of the four bands are presented in Paper I. 

The PBCD fits files are calibrated in terms of MJy sr" 1 
and the photometry tasks used in iraf are designed to 
handle counts/sec. Therefore a proper conversion from the 
flux units to counts and estimating the photon noise is 
required. Photon noise was computed by using a phpadu 
(photons per image count), where phpadu = (exptime x 
gain) / (fluxconv x pixratio). exptime is the total exposure 
time, gain is 3.3, 3.71, 3.8 and 3.8 and fluxconv is 0.1125, 
0.1375, 0.5913 and 0.2008 for channels 1, 2, 3 and 4, re- 
spectively. The pixratio is the ratio between the area of the 
original IRAC pixels and the pixels in the mosaics. This was 
1 in our case. 

Aperture photometry of point sources was extracted us- 
ing an aperture of 2.4". Sky correction was applied using an 
inner sky annulus of 2.4" and an outer sky annulus of 7.3". 
Zero-point magnitudes corresponding to these aperture set- 
tings (including the aperture corrections) were 17.79, 17.30, 
16.71 and 15.89 respectively for Channels 1, 2, 3 and 4. Us- 
ing the individual photometry we then obtained matched 
sources between the various bands. A total of 1580 point 
sources were matched between all four bands. The mean 
photometric uncertainities are 0.02, 0.02, 0.07 and 0.15 mag 
for channels 1, 2, 3 and 4, respectively. 



3 RESULTS 

3.1 Extinction maps 

In Fig. 2, the grey scale background shows the near-infrared 
extinction map of the region. This map was produced by 
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Figure 3. A zoom-in view of the region around DR21(OH) and 
the associated dense filament, with the three colour composite 
background shown in Fig. 1 overlaid by 6th nearest neighbour 
stellar density contours. Note the structure in clustering and the 
stars that form spokes cluster. 



Nyquist binning the spatial distribution of extinction val- 
ues with a bin size of 150". The JHK photometric data 
were plotted on an H-K vs J-H colour-colour diagram. In 
this diagram (not shown), all stars were deredenned to a 
locus perpendicular to the redenning vector and the vector 
length was measured appropriately converted to extinction 
values. Given the large distance to this region, the observed 
area is crowded by field stars redenned only by interstel- 
lar extinction. Including these field stars affects the con- 
trast of the extinction map. Therefore, only those points 
with A „ ^5 mag were chosen to produce the extinction map 
shown in Fig. 2. This map traces regions of extinction be- 
tween A„ of 5-30 mag, limited by the completeness limit 
of the WFCAM observations. However, the DR21/W75 re- 
gion are associated with several dense regions with extinc- 
tion greater than A„ ~30 mag which do not show any stars 
at 2 ^m, resulting in missing data points in our map. These 
pockets appear as white regions adjacent to the black spots 
(for example see the white patch adjacent to the black head 
of the region marked "C" in Fig. 2). Such white patches 
are supplemented appropriately for the portion of the field 
covered by SCUBA observations. This was accomplished by 
converting the SCUBA column densities to extinction values 
and an absolute calibration made by using robust extinction 
data points that was measured with NIR data. In those re- 
gions where SCUBA data is not available the white patches 
remain. A histogram of all the extinction values obtained 
from the WFCAM data suggest that the mean extinction 
to the DR21/W75N regions is A„=6-8 mag. This i s in good 
agreement with the value obtained bv lJoshil j2005l) towards 
this line of sight within the galaxy. 

The extinction map clearly traces the north-south dense 
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Table 1. Embedded Clusters in the mapped region 



Region 


Coords (J2000)(deg) 


Radius (pc) 


Density (stars pc 2 ) 


Isoperimetric 


T 


name 


a 


S 


effective 


core 


half-power 


peak 


quotient 




W75N 


309.652 


42.6462 


0.96 


0.63 


20 


44 


0.76 


0.66 


Diamond ring 


309.660 


42.1197 


1.62 


1.62 


27 


54 


0.20 


0.99 


L906E 


309.383 


42.2238 


1.24 


0.91 


23 


49 


0.45 


0.73 


DR21/DR21(OH) 


309.767 


42.3764 


1.46 


1.08 


24 


48 


0.23 


0.74 



k-33 
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Figure 4. Colour-Colour (CC) diagram from the Spitzer IRAC 
data for DR21/W75, including stars associated with the Diamond 
Ring region and L906E. The solid box represents the region of 
Class II sources (Allen et al 2004). The dotted line shows synthetic 
fluxes for a YSO with a temperature of 32000 K; the solid curve 
shows a YSO track at a temperature of 15000 K (Whitney et al. 
2004). Numbers and letters refer to specific sources discussed in 
Paper I. 

molecular filament encompassing DR21, DR21-OH and the 
W75N region (Paper I). It also reveals other regions of high 
extinction towards the diamond ring cluster, L906E and re- 
gion "C", which is the darkest patch in Fig. |3 extending 
beyond the limits of our map. It can be noted that these 
regions of high extinction coincide well with the brightest 
emission regions in the MIPS 24^im image of Marston et 
al. (2004) and the 850 /im dust continuum cores in Paper I. 
However, the Spitzer/IRAC images and 850 /^m data does 
not encompass region "C" , so the associated star formation 
activity is not discussed here. 



3.2 Embedded Clusters 

The star forming regions DR21 and W75N appear as bright 
near-infrared nebulae owing to illuminated gas and shock 
excited outflows. However, no significant stellar density en- 
hancements are observed at 2/xm (in the K band) in these 
regions. This means that the embedded population must 



be hidden inside the dense molecular cloud at higher ex- 
tinction or that the YSO population is too young to be 
bright at 2/^m. This supposition becomes obvious with the 
availability of IRAC images from the Spitzer Space Tele- 
scope, where many embedded YSO's are clearly visible. The 
contours in Fig[2] trace the embedded clusters in the re- 
gion as identified by stellar density enhancements in the 
4.5/im point sources. These contours were obtained by com- 
puting the 20th neares t neighbour (20NN) en hancements 
llCasertano fc Huf] |l985: Fcrr eira fc Lad a 2006) of the spa- 
tial distribution of 14,063 sources detected in this band. The 
clusters were detected by using the 4er contours above the 
mean star counts in the region where a is the variance of the 
mean star count values. Table. 1 lists the embedded clusters 
identified and their properties. Embedded clusters associ- 
ated with DR21/DR21(OH), W75N, the diamond ring re- 
gion and L906E were detected. There are a couple of fainter 
clusters marked clusl and clus2 in Fig. [5] but these do not 
posess significant peaks, so they are not listed in Table. 1. 
The embedded clusters are found to be spatially coincident 
with the high extinction regions revealed by the extinction 
map obtained from the near-infrared WFCAM data. In fact, 
the observed clusters are better aligned with the whitish 
patches which represent missing data points in the near- 
infrared extinction map. The DR21/DR21(OH) cluster also 
aligns well towards one edge of the dense filament traced by 
the 850 /im emission. 

The 4 a contour was used to identify the boundary of 
the cluster, its effective radius and isoperimetric quotient. 
To compute the peak positions, core and half-power den- 
sities and the r parameter we made use of a 20NN map. 
Table. 1 lists effective radius, peak and average stellar den- 
sities, and the r parameter for the clusters. The effective 
radius is computed by measuring the area A enclosed by the 
4 a contour and computing R—y/A/n. In Figure. 3 we show 
contours of a 6NN map to demonstrate the structure within 
the clusters; it can be noted that the clusters are elongated 
in the direction of the dense filament and they display signif- 
icant sub-structure. Two parameters namely the isoperimet- 
ric quoti ent and r are used to measure the structure of the 
clusters l|Ferreira fc Ladall2006l ) . The isoperimetric quotient 
for each cluster is defined as follows: If the cluster bound- 
ary contour has a perimeter p units enclosing an area A 
square units, then the isoperimetric quotient of this curve 
is 4-irA/p 2 . It measures the ratio of the area enclosed by 
the cluster boundary to the area of a circle with the same 
perimeter. If our clusters were circular, this quotient will 
be equal to 1. r is defined as the ratio of the core radius 
to the effective radius and measures how much of the clus- 
ter population is in th e core and how much is in the halo. 
iFerreira fc Ladal i2006l) defines a cluster as heirarchical if r 
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Figure 5. Colour magnitude diagrams from near-infrared JHK 
data. In both panels, points represent data derived from WFCAM 
photometry. The solid line shows lMyr old PMS model track from 
Baraffe et al., the dotted line represents the 3Myr old ZAMS track 
from Geneva models and the solid arrows represent redenning 
vectors of A„ = 30 mag. 
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Spectral Index (Observed) 

Figure 6. A histogram of spectral indices measured from linear 
least-squared fits to source photometry in the four IRAC bands. 
Thus, only sources detected in all four bands are plotted. A similar 
plot is presented for the IC348 Spitzer observations of Lada et al. 
(2006) for comparision. 



is greater than 0.455. The last two columns in Table. 1 lists 
the isoperimetric quotient and r values for each cluster. It 
can be noted that the DR21 and diamond ring clusters are 
the most elongated and heirarchical, while the W75N cluster 
is the most circular and centrally condensed. 

3.3 Color-Color(CC) and Color-Magnitude(CM) 
diagrams 

The nature of the embedded young stellar population can 
be evaluated using the photometric data plotted on colour- 
colour (CC) and colour-magnitude (CM) diagrams. The ma- 
jority of the embedded population is bright in the IRAC 
bands. Therefore, the 1,592 point sources matched between 
all four IRAC bands were used to plot the IRAC CC and CM 
diagrams. Fig. 2] shows a [5.8-8.0] vs [3.6-4.5] CC diagram. 
The numbers and letters in these plots refer to specific tar- 
gets that are discussed in detail in Paper I: the numbers are 
bright 8.0 /im point sources, often compact clusters and/or 
luminous sources in the region, while the letters are can- 
didate outflow driving sources. Almost all exhibit colours 
consistent with protostars. In Fig.^Jthe redenning vector is 
shown with an arrow, the two cur ves represent Class I mod- 
els for two different temperatures i Whitney et all2004i) . The 
dotted curve represents the colours of a Class I YSO at a 
temperature of 32000 K for various inclination angles, and 
the solid line is a similar object at a temperature of 15000 K. 
The square box represents the zone of Class II objects as 
discussed bv lAllen et al] i2004h . The concentrated circular 
distribution of points to the lower-left of the square box 
represents the pure photospheres. Fig. [I] shows that about 
50% of the point sources detected in the IRAC bands are 



pure photospheres, while the remaining are distributed in 
the zone s representative o f younger objects such as Class 0, 
I and II i Allen et al.ll2004l : iMetreath et alJl2004h . Notice the 
few points with excess red colours that occupy the upper 
right-hand regions of the plot, i.e. regions redder than the 
model curves of 15 K and 32 K Class I source. These are 
representative of candidate massive stars. 

Although the JHK data do not show any significant 
embedded clusters, they do reveal several embedded sources 
and the majority of the halo population distributed outside 
the densest regions of the molecular cloud. Also, the higher 
spatial resolution of the near-infrared data can separate mul- 
tiple components (when these are present) among the em- 
bedded population. As mentioned earlier, the JHK merged 
catalog contains ~43,000 sources and ~60,000 sources in the 
H and K bands alone. These points are plotted on two CM 
di agrams shown in Fig. The pre-main-sequence models 
of iBaraffe et al] lll998l) (solid line) and the zero-age-main- 
sequence model trac ks from the Geneva Observ atory stellar 
models (dotted line) <lLeieune fc SchaererlboOlD are plotted 
for reference. The observed points were compared to model 
data with ages between 1 Myr-3 Myr and distances ranging 
between 2-3 kpc. We find that the near-infrared data are 
matched relatively well for an age of 3 Myr and a distance 
of 3 kpc. Note that the age of 3 Myr applies only to the 
population of stars detected in the NIR JHK bands; it does 
not apply to the deeply embedded population revealed by 
the IRAC data. The YSO population traced by IRAC data 
consists largely of younger objects, with an age of perhaps 
1 Myr or less. 
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3.4 Spectral indicies of the embedded population 

The IRAC bands are well suited to measuring the infrared 
excess emission that can be attributed to circumstellar en- 
velopes and disks. The IRAC band fluxes have been suc- 
cessfully used to determine the disk excess objects and ob- 
tain reliable estimates of the disk fracti on in nearby low 
mass star forming r egions such as IC348 jLada et al]l2006ft 
. lLada et al.l <|2006T) show that the pure photospheres and 
the disk excess sources clearly separate out when the IRAC 
fluxes are used to estimate the power law index of the ob- 
served spectral energy distribution. This is because the disk 
emission dominates in the 4.5-8.0/im, where the SED's will 
significantly deviate from the pure photospheres of low mass 
stars. A simple least squares linear fit to the observed SED 
between 3.6pm and 8.0pm was obtained for all of the 1,580 
sources found in all four IRAC bands. Table. 2 (available 
online) lists the magnitudes in the four bands and the com- 
puted spectral indices. The SED's were not deredenned be- 
fore fitting because an accurate estimate of the A„ to each 
source is possible only by fitting model atmospheres. This 
can not be done with the limited data available for this re- 
gion. 

A histogram of the observed spectral indicies a for all 
sources detected in the four IRAC bands in DR21/W75 is 
shown in Fig. For comparision we also plot the indices 
measured for IRAC sources observed i n the low m a ss sta r 
forming region, IC348, as described bv lLada et, al.l |2006). 
DR21/W75N is located at 3kpc resulting in a large line-of- 
sight extinction that will stretch the a values. Extinction 
due to a fixed column density of gas and dust is a function 
of wavelength which will attenuate the shorter wavelengths 
much more than the longer wavelengths. This will result in 
an artificial rise in the slope of the spectral energy distribu- 
tion if the foreground extinction is significant, resulting in a 
higher value of a. However, the change is larger for sources 
with lower a and relatively smaller for sources with higher a 
thus causing a stretching of the a values rather than a sim- 
ple increase. We therefore corrected the IRAC ma gnitudes 
for su ch effects by using the standard extinction law iMathisI 
1990) and assuming an uniform foreground extinction value 
of Av=9mag. Note, however, that it does not correct the ac- 
tual extinction to each source which is due to the molecular 
cores or disks. In DR21/W75, the observed spectral indicies 
range from -4.0 to +4.3. The histogram for IC348 shows a 
prominent peak at a ~-2.7, which corresponds to pure pho- 
tospheres, and a somewhat less prominen t peak at a ~-1.8 , 
corresponding to sources with thick disks jLada et all200fj) . 
Similar peaks can be seen for the DR21 data as well, the 
shapes of both histograms compare quite well, representing 
similar aspects of star formation in both regions. However, 
the histogram for DR21 deviates from that of IC348 at a 
values greater than and particularly beyond 2. A small, 
but significant fraction of sources in DR21 are found to have 
high spectral indicies, between 2 and 4.5. The higher values 
of a is not merely an effect of the stretching effect explained 
above due to high density regions. The high a sources are 
better represent the source functions because the change 5a 
is smaller for high a sources compared to low a sources. 
We shall discuss the implications of this result in the next 
section. 

In Fig. |7|we show the spatial distribution of the IRAC 



point sources using circular symbols. Sources with a >3 are 
shown by cross marks. The size of each symbol is propor- 
tional to the observed spectral index of the corresponding 
source. The background greyscale image is a 4.5pm mosaic 
image of the region. Sources with the highest spectral index 
are concentrated in the densest regions associated with the 
north-south molecular ridge that runs through DR21 (see 
Fig. 1), W75N, L906E in the west-south- west, and the Dia- 
mond Ring region. In particular, strong clustering of sources 
can be noticed around DR21, DR210H, W75N and L906E 
region (see Fig. and Fig. 7). These clusters are connected 
by a distributed population of sources with high spectral in- 
dices, that can be found spread along two filaments, one run- 
ning north-south and the other towards the south-east and 
north-west. A general correlation of the disk sources with 
the dark patches in the entire field can be seen in Fig. [7] 



4 DISCUSSION 

We have explored the stellar popula tion in and a round 
two major, lumino us (~ 1O 5 L0 iMoore et alJ Il988l : 
ICampbell et alJll982>) star forming regions, namely DR21 
and W75N, using wide field UKIRT/ WFCAM JHK and 
Spitzer/IRAC 3.6, 4.5, 5.8 and 8.0pm images. The ob- 
served area encompasses several l uminous far-infrared 
sources, clusters of UC HH regions jTorrelles et alJ Il997l : 
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tion. The estimated distance to all these sources is more 
or less the same; the analysis of the data presented in 
Sec. 3.3 suggests a distance of ~3 kpc, which is in good 
agreement with other estimates derived from various 
methods fCenzel fc Downesl Il977t ICampbell et al.l Il982l : 



iFischer et alJll985UOdenwald fc Schwartzlll993F 

The results presented in the previous section demon- 
strate that the DR21, DR21(OH) and W75N regions contain 
young stellar objects that can be classified as Class 0/1 or 
Class II, based on Spitzer colours and spectral indices. This 
implies that star formation is a recent phenomenon in all of 
these regions. The majority of the young stellar population 
in each area is so deeply embedded that it is not uncov- 
ered in observations at wavelengths shortward of 2pm. The 
SCUBA 850 pm cores described in Paper I and shown here 
in Fig. 1 probably represent an evolutionary phase that is 
earlier even than that traced here with Spitzer. The SCUBA 
cores (32 are identified in Paper I) are typically not associ- 
ated with 8.0 /im sources. They are generally confined to a 
narrow, north-south filament in DR21, and to a compact 
region around W75N. The Spitzer YSOs, though also cer- 
tainly clustered around the dense cores and filaments, are 
none-the-less somewhat more widely distributed. This could 
be due to the displacement of YSOs away from the central 
density peaks as they evolve. Even so, the different epochs 
of star formation traced by WFCAM, Spitzer and SCUBA 
all seem to be co-existing in DR21/W75. Moreover, com- 
mon trends such as distance, YSO class and estimated age 
suggest that the individual sites of star formation - W75N, 
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Figure 7. Distribution of infrared excess sources. IRAC channel 2 image is shown in grey scale and the circular symbols represent the 
excess sources. Crosses represent sources with a >3. The size of the symbols are proportional to the value of the spectral index. 
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DR21, DR21(OH), the Diamond Ring and L906E - may be 
connected in a global way. In the following we shall argue 
that there is more evidence to support this scenario. 

The 8^tm channel of the IRAC instrument traces warm 
dust and emission from PAHs (polycyclic aromatic hydro- 
carbons) much more efficiently than the shorter wavelength 
bands. The colour image of the entire region shown in Fig. 1 
displays prominent filamentary nebulae. At least five long fil- 
aments of nebulosity (marked by the sketched lines in Fig. 1) 
can be identified in the observed region, four of which con- 
verge on the DR21/DR210H regions. Note that the bases of 
these filaments are traced to some extent even in the 850 /im 
emission, showing that the filaments are indeed made of 
dense gas. We have seen from Fig. |7| that a visibly promi- 
nent distribution of thick disk sources is concentrated along 
the thickest filament coinciding with the DR21, DR2IOH 
and FIR 1/2/3 regions, although the W75N region is well 
away from these filaments and appears disconnected. A care- 
ful examination of the same figure will also show that two 
other filaments that originate at DR21 and run towards the 
northwest and a third directly westward of DR21 are also 
correlated with the distribution of thick disk sources. The 
thin filament running directly westward of DR21 and linked 
to the outflow may be due to the outflow punching out of the 
filament. While the correlation of embedded YSOs and fil- 
aments is prominent for the filaments converging on DR21, 
the same is not true for the diamond ring region. The dia- 
mond ring region has a bright filament running northwards 
that curves toward the L906E region marked in Fig. 1. The 
YSOs in Fig. 7 all appear to be concentrated in the L906E 
region and just below the filament, rather than at the head 
of the diamond ring. This suggests that the diamond ring 
region is relatively evolved, while L906E is actively form- 
ing young stars. Note that no H2 jets were discovered near 
the diamond ring, while a number of flows were identified 
around L906E (Paper I). A proper analysis of the filaments 
and their relative orientation requires additional mapping 
at 850 fj,m (our SCUBA map does not extend as far south 
as the diamond ring, or as far west as L906E) and CO line 
data with good spatial and velocity resolution. 

As shown in section. 3.2, embedded clusters are asso- 
ciated with each of the main targets in the region, namely 
DR21/DR21(OH), W75N, L906E and the diamond ring re- 
gion. The clusters show cores on scales of ~1 pc. If each clus- 
ter peak represents a molecular clump out of which the clus- 
ters were born, then the size scales of the cluster cores rep- 
resent the scales of fragmentation in this filamentary cloud. 
The main DR21 filament traced by the 850 /im emission is 
about 16 pc long and is one of the most structured fila- 
ments in terms of embedded clusters. A zoom-in view of 
this region is shown in Fig. 3 along with contours of the 
stellar density. The contours display multiple stellar density 
peaks. These clusters are elongated roughly along the north- 
south filament. In this figure notice the bright stars marked 
as "Spokes Clusters", which include the FIR 1/2/3 850 /im 
peaks (refer Fig. 1) in one line. Each of these bright stars 
are surrounded by a noticeable linear alignment of fainter 
stars that are arranged like the spokes of a wheel whose 
centre coincides with the bright sta r. This type of struc- 
ture has been recently identified by iTeixeira et alJ (I2006T) 
in NGC2264, which is at a somewhat closer distance of 
d=800 pc, and is therefore more clearly visible. These au- 



thors show that the linearly aligned fainter stars are indeed 
redder than the central luminous source. Such structures are 
thought to be important because they represent primordial 
structures in the f ormation of a star cluster representin g 
true proto-clusters llKurosawa et al 1l2004l : iBate et aljEo 03 ) . 
Although the spokes configuration in the clu sters shown in 
Fig. 3 are not as prominent as that shown bv ITeixeira et all 
(2006) (possibly due to the larger distance of DR21), there 
are never-the-less four such clusters in DR21. The associ- 
ation of these clusters with bright 850 /im emission and a 
young population of stars with high spectral indicies indi- 
cates that this region may represent multiple adjacently lo- 
cated protoclusters. 

The filamentary structures found in the DR21/W75 re- 
gion on the scales of giant molecular clouds (16-20 pc), and 
those found on smaller scales, such as the "spokes clusters", 
collectively resemb le the patte r n pre dicted by the numer- 
ical simulations of IBate et al.l (l2003f) . The simulations by 
IBate et alJ J2003ft are made for a 5OM0 cloud, which is much 
smaller than the actual mass of the DR21/W75 clouds. How- 
ever, qualitatively, their results seem to be applicable even 
here, suggesting that the mechanism of gravo-turbulent star 
formation effectively operates on the scales of giant molec- 
ular clouds, up to several hundred thousand solar masses 
and several parsecs in dimension. The filamentary nature 
appears heirarchical if we consider the smaller scale spokes 
configurations and the larger scale filaments feeding on to 
the center of the DR21/DR21(OH) region. This heirarchy 
can be carefully observed by comparing structures in Fig. 1 
and Fig. 3. 

Finally, we mention the implications of the sources with 
high a values in DR21 distributed along the dense filament 
and also coinciding with sign-posts of massive star forma- 
tion, a >2 in a simplistic interpretation implies p resence of 
Class objects. Recently. iRobitaille et all i2006t) presented 
and discussed a large grid of 200,000 YSO model SEDs. Ac- 
cording to their work, there are a number of reasons why 
a (measured in the Spitzer IRAC bands) can assume high 
values: a) younger sources such as Class sources are known 
to have high a, b) higher mass accretion rate and/or larger 
disk/envelope masses boost a, and c) beyond 5000K, a will 
increase with the temperature of the star, due to a lower 
contribution from the stellar flux and a higher contribution 
from the infrared dust spectrum. In DR21/W75N the ob- 
served data points satisfy almost all of the above conditions 
because they are situated in the densest parts of the star 
forming region, they coincide with sign-posts of massive star 
formation such as UCHII region and water masers, and be- 
cause they occupy the colour zone of Class protostars on a 
CC diagram. More quantitative comparisions of the real na- 
ture of massive stars requires the plotting of data from theo- 
retical SEDs on to CM diagrams. However, in colour-colour 
space there is significant overlap of data points for the differ- 
ent evolutionary states, sour ce masses, temperatur es, disks 
and even inclination ang les (IRobitaille et al.l l2006). partic- 
ularly if one uses a limited number of photometric bands. 
Rigorous comparisions will require multiple wavelength ob- 
servations, particularly at larger wavelengths such as 24 /im 
and 100 fj,m where younger and luminous objects are known 
to particularly bright. 
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5 CONCLUSIONS 

We analyse NIR and MIR point-source photometry over an 
area that includes the massive star forming regions DR21, 
DR21(OH), and W75N, as well as recently-discovered clus- 
ters of more evolved young stars (the Diamond Ring re- 
gion), and bok-globule-like cores where low mass stars are 
probably forming (L906E) . We use colour-colour and colour- 
magnitude diagrams, as well as extinction maps and sur- 
face density plots, to examine the population of young stars 
across the region. We also measure the MIR spectral index 
for 1580 sources. 

Our extinction map traces much of the clumpy, north- 
south filament associated with DR21 and DR21(OH), as well 
as new cores that may be associated with separate episodes 
of star formation. The very densest cores, however, are not 
extracted from the NIR data, because of extinction larger 
than A v ~30 mag. 

NIR and MIR color-colour and colour-magnitude plots 
show that the JHK and 4.5 pm-8.0 fim data trace differ- 
ent YSO populations, as expected: the WFCAM data re- 
veal sources with ages of ~l-3 Myr, and mostly present in a 
widespread halo, while the IRAC point sources are clustered 
along the dense filaments traced by 850 /jm data. Spectral 
indices for the Spitzer sources indicate a more embedded 
and/or massive population in DR21/W75 than is seen in 
the IC348 low-mass star forming region. There is a signifi- 
cant fraction of sources with indices greater than 2 and up 
to 4 in DR21/W75. Indeed, ~40% of the Spitzer-detected 
sources in the field have SEDs consistent with Class 0/1 or 
Class II YSOs. The reddest sources are generally associated 
with the densest cores, as traced in 850 /im dust continuum 
emission in Paper I, or in 8.0 (jm PAH emission. The NIR 
data, tracing the evolved and halo YSOs represent a ~ 3 Myr 
old population whereas the Spitzer data, tracing the more 
deeply embedded content and also associated with outflow 
activity, represents a relatively younger population of 1 Myr 
or less. 

Star formation in DR21 region is found to be occuring 
in filaments of 10-20 pc lengths that is fragmented in to mul- 
tiple clusters on the scales of 1-2 pc. We identify four spokes 
cluster configurations in the main DR21 filament that is rep- 
resentative of proto-clusters. The morphological appearance 
on the scales of the entire Giant Molecular Cloud resem- 
bles well the num erical simulations of cluster formation by 
iBate et all (|2003h . 
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